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Abstract 



The paths of stars in galaxies have circular velocity independent of their distance 
from the centre of the galaxy. Newtonian mechanics with a logarithmic potential has 
such paths. In relativity these paths can be taken to be geodesies and this requirement 
places a resticting equation on the spacetime metric. This equation has a non-unique 
■^j- , spherically symmetric solution that in the newtonian limit has a logarithmic potential. 

ON ' It also can be solved in terms of a conformal factor. In particular it has solutions which 

are conformal to the vacuum-Einstein solutions and such spacetimes are solutions to 
the vacuum-Bach equations. Therefore it is suggested that the Bach equations describe 
dynamics on galactic lenght scales. 

q 1 Introduction 

The majority of galaxies conform to a few stereotyped structures such as spirals and barred 
spirals. The speed at which stars rotate around the centre of such galaxies is not what 
would expected from their detected mass - rather the rotation speed is independent of the 
distance of the star from the centre of the galaxy. Both optical and radio observations 
support this conclusion. The usual solution to this problem is to invoke the existence of a 
large amount of as yet undetected material which resides in contrived places. This is the 
so-called the "missing mass" problem. Modern relativistic theories of gravitation equate 
the effects of matter with curvature. Curvature is described by the Riemann tensor which 
contains derivatives and products of the Christoffel symbol. The Christoffel symbol is 
constructed from derivatives of the metric, which is an indicator of the distance between 
points. In §2 it is shown that in Newtonian kinematics, constant orbital velocity requires 
a potential of ln(r); this result was mentioned without derivation in ||. The logarthimic 
potential is related to the metric and gives Christoffel symbols with terms proportional to 
1/r. The products and derivatives of these are both proportional to 1/r 2 ; the standard 
approximation scheme uses only the highest derivative terms and hence production of the 
linearized Riemann tensor breaks down. Thus it is impossible to establish contact between 
linearized relativistic theory and Newtonian theory solely from knowledge of a logarithmic 
potential. The relativistic treatment has to be approached in a different manner: the 
method that I suggest is to note that the paths of particles are given by geodesies which 
are curves of minimum length. For galaxies the geodesies are known - they are simply the 
circulating stars' rotation curves, called galactic rotation curves. These have recently been 
reviewed in detail by Blattaner and Florido ||] (2000). Once the geodesies are known it 
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is possible to find a coordinate dependent constraint on the metric of the spacetime that 
the geodesies inhabit; how this is done is described in §3. 

Observations of interaction between galaxies suggest that even more missing mass is 
required, and in cosmology large amounts of missing mass are also postulated: the pattern 
is the longer the length scale the higher the proportion of hypothetical undetected matter. 
There is a critical density above which the Universe is closed and will contract. Nearly 
all cosmologists feel that the Universe is close to this critical density. Using Einstein's 
equations this critical density requires that about 90% of the Universe is as yet undetected. 
An alternative is that different equations might be appropriate, again the new equations 
originating from a relativistic theory other than Einstein's. It would be expected that 
the appropriate equations will be the same ones that have a solution with a metric which 
explains constant galactic rotation curves. 

There are two main sorts of explanation for constant rotation curves. Currently the 
most widely accepted is that there is a large amount of as yet undetected mass, see 
Blattner and Florido (2001) and references therein. The second is that force laws 
other than Newton's are needed, see in addition to §, Roberts @ (1985) and @ (1991) 
and references therein. There are also explanations built on analogies with other systems, 
such as hurricanes, see Emanuel § (1991). The explanation suggested here is that galactic 
spacetime is a solution to the Bach Q (1921) quadratic equations. Quadratic lagrangians 
can be expressed as two independent parts: a R 2 short-ranged part and a conformally 
invariant C 2 long-ranged part, so that for present purposes only the C 2 part is of interest. 
Specifically varing the C 2 action one obtains the Bach tensor 

B a b = ^C a cd b R cd + AC a j£ cd . (1) 

Any spacetime with vanishing Weyl tensor C a \, C( i = 0, such as R a b = Ag a b, is a solution 
to the vacuum-Bach equations |l|; in particular any conformally flat spacetime and also 
Kottler's solution 

, <? ( _ 2m Ar 2 \ l9 / 2m Ar 2 \^ l9 9 , „ 
ds = ~ 1 - — + — ) dt + 1 - — + — dr + r dS 2> ( 2 ) 



are solutions. This suggests that the large-scale vacuum, see the conclusion of [ 10 1 , can be 
thought of as having a conformally flat metric rather than a flat one. One is left with the 
situation that large scale dynamics are described by the Einstein-Bach equations, the Ein- 
stein part dominating on lenght scales of the solar system and the Bach part dominating 
on galactic lenght scales. This suggests that on cosmological scales field equations corre- 
sponding to quadratic lagrangians provide a better description than using the Einstein field 
equations: cosmology and quadratic lagrangians have been reviewed in Querella || (1999) 
and references therein. How constant rotation curves constrain axi-symmetric spacetimes 
has been discussed in Roberts 1| (1986). What happens as r — ► oo is a problem in general 



[11 1 and particularly with any explanation of constant rotation curves. A solution is to 
take the spacetime solution to be the correct model for what was previously taken to be 
the galactic halo, and this should be used for models of gravitational lensing. If one takes 
the Bach equations to be the correct ones, one simply assumes that for sufficiently large r 
a different solution applies. The situation is similar to stars where one can take a different 
Schwarzschild exterior solution for each star, for galaxies one takes a different anti-deSitter 
solution for each galaxy; a difference with this analogy is it is possible to find the size of 
the coupling constant by linearization for the Schwarzschild solution, but unfortunately 
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for galaxies there is as yet no known linearization to fix any of the quadratic coupling 
constants. Constant rotation curves by themselves do not explain the Tully-Fischer [12] 



(1977) luminosity to speed relationship, but perhaps there is a metric which explains both. 
It is surprising that spherically symmetric spacetimes descibe constant rotation curves at 
all, galaxies rotate and so are better modeled by axi-symmetric spacetimes. 

The conventions used are those of |J|, in particular the signature is (—,+,+,+). Here 
all work is done in the plane, so that 6 = 0, the reason that this can be done without 
loss of generality for spherically symmetric spacetimes are the same as those given in 
Chandrasekhar @ (1983). 

2 Newtonian Analysis. 

In this section it is shown that in Newtonian theory a logarithmic potential gives circular 
paths independent of r. In spherical coordinates Newton's second law is 

Va/m = (V r , V+ym = (r - rift 2 , ^r 2 0), (3) 

by spherical symmetry the (ft component integrates immediately to give 

r 2 (ft = L. (4) 
Subsituting using u = 1/r and using u' = du/dcft and ||, gives 

f = -Lu, r = -L 2 u 2 u". (5) 

Using [5| the velocity vector is 



L(-u',l), (6) 



and also using || the first component of Newtons second law || in polar coordinates is 

d 2 u | u= V r 
deft 2 mL 2 u 2 

which is called the Binet equation: an imporatnt thing to note is that there is an undif- 
ferentiated u on the LHS. Recalling gM> = 1/r 2 ; the square of the velocity at any point 
is 

v 2 = r 2 + (r^) 2 =L 2 (u l2 +u 2 ), (8) 
for circular orbits v! = so that this is 

^ 2 =^ 2 = ^> ( 9 ) 



for circular orbits u" = so that from [7| 



Vr r 3 



mW m 

substituting in |9| 



L2 = -373 = --^n (10) 



1? r 

V 2 c = — = V r . (11) 
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For constant circular orbits 

mv 2 c = -rV r , (12) 

or 

V = — mv 2 ln(r), (13) 
the desired result. The velocity and force vector fields are 

v = (0,L) = (0,rv c ), Yl=(-^,o). (14) 

m \ r J 

To see that the correct rotational velocity /r 2 , rather than just i>^; note that 

using v 2 ^ one gets 

mv 2 , 

4 = L ^ V = -^, (15) 

for || and ^ respectively, the potential is shorter ranged potential than that of Newtonian 
gravitation. 

3 Spherically Symmetric Geodesies. 

Here geodesies are derived by the Lagrange method, as used for example in Chandrasekhar 
H (1983). The lagrangian for geodesies is 

2C = g u i 2 + g H ft + g„r 2 + g ee 2 = p a p a , (16) 

in the (—,+,+,+) convention, 2C = —1 for timelike geodesies. The Euler equations are 

dp a _ dL 

dr dx a ' { ' 



Fori! 

r> ~ = 

dx 



dL ■ ■ ■ ■ • 

p<* = ~n 7 ^ = (9ttt,g<i>ci><t>,grrf,gee0), v a = (t,<t>,r,0), p a p a = -l, (18) 



and 

dC 

2 dx^ = ^ °' 9tt ^ 2 + + 9rr ^ + 3ee,r0 2 , 0). (19) 

The first two components of the Euler equation integrate to give momentum and its 
derivatives 

Pa = (E, L, g rr r, g ee 9), (20) 
where E and L are constants for each geodesic. For the second two components 

dpr .. . . 2 d Pe n d ■ 

— — = g rr r + g rrr r , — — = = —geev, (21) 
ar ar ar 



working in the plane 6 = so the last of these is satisfied, from 19 and ^ the r component 
of the Euler equation gives 

= 2g rr r + g rr , r f 2 - g t t, r i 2 - g^,A 2 , (22) 
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in terms of u" , r takes exactly the same form as in the Newtonian does <b in 

„2 



equation [|. Substituting for t using 21 then dividing ^ by —2Lu and taking 
this can be seen to be a generalization of |?] the Newtonian Binet equation 



— QrrU Qrr,i 



u 12 g tt ,r ( 



E v 



2u 2 



2 V^5tt 



(23) 



The dynamics of geodesies in the plane are described by equations [R] and |22|. For 
circular orbits again u" = u' = =>■ r = r = 0, so that the two equations become 



= g tt J 2 + 



y\ o = 5 «i 2 



2£. 



For non-null geodesies 2C ^ 0, then these equations El have solution 



■2 2£ 

4 ' = T 



2C 

—r-9u,r, h = 9ttg<t><t>,r - 94>4>9tt,r T 0, 



(24) 



(25) 



for null geodesies 24 implies h = 0, so that these equations do not hold. Defining the 
angular velocity as 

It (j) g tt L 
E' 



n = — = 

dt dt/dr t 
so that the geodesic momentum vector ^ is 



(26) 



(27) 



From [25| and |26j the differential equation governing angular motion is 

gtt,r + g<jxf>,rtt 2 = 0. (28) 
For axi-symmetric spacetimes a much longer calculation shows that this generalizes to 



^ 9tt,r [gtt,r — ^9t4>,r + ^ 9<f><p,r) \9tt,r + ^9t<t>,r + ^' 



0. 



For constant angular velocity, integrating 28 

g tt = -ti 2 J g^ r dr = k-VL 2 
taking the constant of integration k = — 1 and g^ = r 2 



gtt 



■(i + flV). 



(29) 



(30) 



(31) 



3 7 J 5 



Note that for the signature (+, - 



remains the same, and the constant of inte- 
gration is k = +1- 31 is of the same form as the de-Sitter metric [2] with A = +3f2 2 . For 
constant rotational velocity [| becomes v 2 = r 2 Q 2 so that 



r gtt,r + 



r v 2 c = 0, 



(32) 



Integrating with = r 2 and, as before, choosing that constant of integration to give the 
correct flat spacetime limit gives 



ds 2 = - (l + 2t> 2 ln(r)) dt 2 + g„dr 2 + r 2 cE 



2 

2) 



(33) 
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which agrees with the newtonian potential [l3|. The tensors for logarithmic metrics of this 
form suggest no easy exact solution. 



This is not the unique solution of 32 by any means, assuming a conformally flat form 
for the metric 

ds 2 = E 2 ds 2 M = E 2 (-dt 2 + dr 2 + r 2 dY>fj (34) 



and after dividing by 2rE, |32| becomes 

(v 2 - l)rE r + v 2 E = 0, (35) 

which has solution 

E = Ar^ (36) 
where A is a constant. Defining the luminosity coordinate 

R = Er = Ar 1 ^ (37) 

the metric becomes 

ds 2 = - A ^-vl)tfvl dt 2 + v 2 c) 2 dR 2 + R 2 d ^ (3g) 

as r — > oo, gu — ► oo so that the metric is not asympotically flat. The metric ^ is 
conformally flat so that it is a solution to the vacuum Bach equations. A generalization 
of |34| is to metrics conformal to known solutions, c.f.eq.2.28|5|. The generalization of |35| 

Er ^ 9tt ' r (39) 



E (v 2 + 9tt)r 2{v 2 c +g tt y 
which has particular solutions 

.2 

1/3 



Ar^{-v 2 -g tt } (40) 



where = l/(2(v 2 - 1)) for deSitter spacetime and (3 = (3v 2 - 1)/(2(1 - v 2 )) for 
Schwarzschild spacetime and there is no for both m, A ^ 0. The resulting line ele- 
ments are not explicitly calculable in terms of the luminosity coordinate R. 

Another approach to constant rotation curves is to assume a spherically symetric solu- 
tion and use p9 to turn the spacetime an axially symmetric one with the correct geodesies. 



The first bracketed term of 29 integrates to 



If m Ar 3 
gt<p = v c r + — / g t t, r rdr = v c r ln(r) + - — , (41) 



[label g24] where the integral is evaluated for Kottler's spacetime ||. 

4 Conclusion. 

The problem of constant galactic rotation curves has been looked at in three ways. The 
first suggested metric |33| agrees with the newtonian solution to the problem, however its 
tensors appear to obey no simple field equations, as there remains an arbitrary function 
g rr in the line element this might be rectified in the future. The second is to produce 



spacetimes conformal |3_| to known solutions, these spacetimes obey the vacuum-Bach 
equations |l[ but in approximation do not give back the newtonian theory. The third 
is to add an axi-symmetric component to known spherically symmetric solutions, these 
spacetimes do not seem to obey known field equations or reduce to newtonian theory. 
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